N 75 20385 

NOTE: This paper is as it was presented during the workshop. 

The Author's review and editorial comments were not re- 
ceived. His slides and figures appear at the end of this 

session. 

THE AEROTHERMAL ENVIRONMENT AND MATERIAL RESPONSE, A REVIEW 

William E. Nicolet 
Aerotherm Acurex Corporation 

MR. WILLIAM NICOLET: Thank you. In response to the letter 

of invitation, as I recall, the wording was that we were invited to 
review and assess current states of technology and make recommen- 
dations, and so I addressed myself to that, rather than giving a 
lot of numbers. It seems like I've been promised, repeatedly, to 
give numbers. I did give a few just to orient tne audience, 
but this will not be a presentation oriented to that end. 

In addition, in the initial response to the letter, I pro- 
mised to review both aerothermal environments and material re- • 
sponse. After looking at the time allocation, I decided I'd better 
delete material response and leave that to this afternoon's ses- 
sions and to other people. So, the focus of this particular talk 
will be a review of the aerothermal environment. 

Figure 5-43 - I'm going to end up duplicating some of the 
material that Jerry presented, clearly, but let's start off by 
looking at the flow and the material response as Aerotherm sees it 
as opposed to how Langley sees it. There are, pretty clearly, a 
lot of overlaps here. 

To begin with, you have a normal shock wave with some re- 
laxation zone behind it, usually of some maximum thickness, at 
the stagnation point. Typically, there is the hot 
shock layer of gases behind it emitting radiation to the body. There 
is some type of a mixing region, hopefully out in the middle of the 
shock layer, bounded by ablation gases flowing inviscidly out from 
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the body on the inside, and by environmental gases flowing more 
or less tangentially on the outside. One might expect this to be 
laminar in the region near the stagnation point with transition 
of turbulence back further. There will be absorption of radi- 
ation in the mixing layer and in the ablation layer. In addition, 
there appears to be important radiation components emitted by the 
mixing layer itself. If one goes over and looks at the other 

end - and I am just going to touch on this - as I said, impor- 

tant absorption in the ablation layer? important events going 
on at the surface: thermochemical events, mechanical removal, 

radiation emission, reflection, melting, depending on the type of 
ablator selected; important events going on in depth: heat con- 

duction, pyrolysis gas formation, scattering; again, depending 
on the material selected. 

Figures 5-44, 45, and 46 are three slides that I will put 
up here really just to allow us to focus down to some numbers. 

To begin with, note that the solid lines are for Saturn, and the 

dashed lines are for Uranus. This is the stagnation-point ra- 
diative heating flux as a function of time. 

To begin with, two different atmospheres are considered here; 

the cold dense and warm atmospheres for both planets. Also two 

different body shapes were considered. Most of the data is for 

a 60° aft angle cone, but the very high radiative flux (above 60 kW/ 
2 

cm ) was computed for an Apollo-type configuration. The convec- 
tive fluxes show slight quantitative differences but, qualitatively, 
are very similar. In contrast, the radiative fluxes are vastly 
different, with the Uranus cold-dense fluxes being nearly an order 
of magnitude greater than those for the Saturn cold-dense entries . 
Moreover, entries into cold-dense atmospheres have radiative flux 
levels which are at least an order-of-magnitude greater than the 
corresponding entries into the nominal or warm entries for the same 
planet. This point will be made over and over again, but has to 
do with the composition of the atmospheres and almost nothing else. 
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RADI A l IVE HEAT BLOCKAGE CORRELATIONS FOR 
CARBON AND Si0 2 ABLATION 
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Figure 5-45 presents the spectral distributions of the inci- 
dent radiation flux. Typically, these are calculated from our 
computer program. These would be for the Saturn nominal entry, 
a relatively steep entry into the Saturn nominal atmosphere. The 
plot on the left shows the radiation coming from the shock layer 
alone . 

On the right we see the radiation with ablation products. 

Note that we have cut out a good part of the radiation in the 
U.V. Simultaneously, we added radiation in the visible portion 
of the spectrum. This is the radiation coming directly from the 
mixing layer. In both of these figures, the clear parts refer to 
continuum radiation; the slashed parts refer to line radiation. 

One might hope to get from a detailed calculation of the ra- 
diation heat transfer correction and blockage correlation of the 
nature shown in Figure 5-46. The solid line represents work that 
was done in support of a Jupiter entry study. It was done three 
or four years ago by Ken Wilson, and subsequently correlated by 
Bill Page. The focus there was for large blowing rates. The 
Jupiter entry case was very severe. Typically, we would see it 
reduce the radiation flux by about a factor of two. My point in 
doing additional calculations for smaller blowing rates which is 
important in the Saturn-Uranus nominal type entries was to inves- 
tigate the effects due to the mixing layer radiation that I dis- 
cussed in the previous slide. As you see, typically, we have im- 
portant additive effects. 

These types of correlations developed from stagnation-point 
solutions, are generally used for the whole body. The objection 
Mr. Walberg was making a few minutes ago was that, in fact, these 
might change shape as you go around the body. John Howell and C. H. 
Liu, here at Ames, are greatly expanding the matrix of calculations 
on this particular subject. Again, it is focused primarily on the 
stagnation point, but it, supposedly, would do a lot to firm confi- 
dence in this type of calculation and the correlations of it. 
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Hopefully, I have set the stage for the type of numbers, 
the type of effects, the type of events we are dealing with, I 
would like to now review calculational and experimental approaches 
to solve the problem, focusing on areas where there are uncertain- 
ties and suggesting, as a last item, ways to reduce the uncertain- 
ties . 


The most uncertain item - and notice on Figure 5-47 I have 
listed input conditions now - the most uncertain item in the whole 
analysis is the atmospheric composition and, particular, elemental 
mass fractions of helium. My calculations jump up and down and go 
all over the place, depending on what we assume there. In par- 
ticular for the Uranus case we can find fantastic radiation fluxes 
for a high-helium-content mixture; and almost none for a low- 
helium-content mixture. 

The elemental mass fraction of the primary radiating species 
in the environment is important, as are the atmospheric scale 
heights. But they certainly take a distant, second place in impor- 
tance to the elemental mass fraction of helium. 

I am going to briefly run down some of the calculational 
methods. To begin with, all of my focus will be on the 2-D flow 
capabilities. There is a figure in the handout dealing with 3-D 
capabilities but, for those of you who are interested in the angle 
of attack, I suggest you look at that and perhaps, talk to me. I 
will not discuss it as part of the oral presentation. 

Let's start with the inviscid type of calculations (Figure 
5-48) applicable right behind the shock front. There are a number 
of finite - difference, time dependent or integral relation methods 
- Jerry Walberg alluded to some - focusing primarily on situations 
where there is no radiation coupling. These would be used for 
basic studies or pressure or boundary layer edge velocities, and 
the like. They would be applicable in the cases where the ra- 
diation is not important. If we add radiation coupling on the 
second line (indicating.) ;, we find that there are a couple of cal- 
culations that can be done. There are a couple of codes available 
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REVIEW OF 2-D FLOW CAPABILITIES 
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with selected organizations. If we go on down to talk about what 
would be required to support heat shield sizing studies, correla- 
tions of shock shape, pressure distributions, and the like are 
certainly vital inputs. They are generally available for the 
shapes of interest for the non-radiation coupling situation. 

Let's go on to boundary or merged layers. Figure 5-49. Here 
my terminology for merged layer is the same that Jerry was using 
for viscous shock layer, that is, a boundary-layer-type calculation 
extending from the shock wave all the way to the body. I tend to 
use them interchangeably since the mathematics tends to be quite 
similar. 

If we talk about a finite difference method coupled to abla- 
tion chemistry, with the laminar or turbulent flow, but without 
radiation, we have such codes as the BLIMP that has been discussed 
previously. It is operational without radiation coupling. There 
are other codes like it, provided that the blowing rates remain 
modest. If we add radiation coupling; same types of codes, same 
types of restrictions. If we reduce, or subtract off, the tur- 
bulent flow requirement we have codes that are applicable for all 
blowing conditions, and this is certainly the situation for the 
typical outer planetary entry of interest. 

If we go on down and ask about a finite difference approach 
considering finite rate chemistry, even without radiation coupling 
and without turbulence, we find that this type of approach has 
generally not been used in the planetary entry situation, although 
the RV community has developed that type of code and some capa- 
bility does exist. I point out that this type of discussion was 
made before I was aware of the most recent presentation of the 
people from JPL. 

Figure 5-50 continues and gets more into intermediate or 
tool-type things that would support heat shield sizing, there are 
various stream tube methods. I would consider Olstad's method a 
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stream tube method. Typically, they can't handle radiation. 

They will do things like radiation cooling factors or cold wall 
heating distributions. If you add ablation they also have some 
capability, but it is more limited. 

If we go into straight integral methods, we can handle 
most of the items of interest, provided we restrain ourselves to 
cold wall events - no blowing or whatever - and we have to resort 
to other types of approximate methods to get the radiation fluxes. 
These are the methods that I have typically used in support of 
heat shield sizing. 

Finally, we have correlations which, again, will be required 
for the heat shield sizing. That would include blowing corrections 
for the convection and the radiation and would refer back to the 
figure I showed before. Some are available for such ablation 
species as carbon and Si0 2 and there are efforts underway right 
now to expand the correlation base. 

I would like to go now to Figure 5-51, review of transport 
properties. This is with application to input to the flowfield 
calculations. To begin, there is a total properties approach, 
and this is a classic approach that has been used for years. It 
was originated by Butler and Brokaw. The entry calculations that 
have been done with it are almost without number. It is very sim- 
ple, however, it is restricted to non-varying elemental composition 
across the layer. And that, in effect, restricts the calculations 
to no ablation or to ablation of a gas which has the same elemen- 
tal composition as the environmental gas. So, with that restric- 
tion, that approach is losing favor. 

There is a series here of three successively more compli- 
cated approaches, namely: correlations for such properties as 
viscosity, diffusivity, thermoconductivity plus equal diffusion; 
coefficient approximation, bifurcation approximation, actual solu- 
tion to the first order Chapman; Enskog solutions. These successively 
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increase the accuracy of the solutions to the diffusion equa- 
tions. They all run into trouble when ionization begins to be 
important. This would, typically, be in the 8,000° to 10,000° 
Kelvin range. It is my feeling that higher order solutions or 
improved correlations are required to go above that substantially. 

Figure 5-52 reviews the radiation transport codes or pro- 
perties that are available, detailed codes have been generated and 
are available from several organizations. They are used in support 
of basic studies, reduction of experimental data. There is at 
least one that is available that will support flow coupling cal- 
culations. Typically, they are also used to define multi-group 
radiation models. 

Concerning the properties, these have been the subject of 
a recent review at Langley. It is my feeling - although I haven't 
read the report yet - that the environmental gases are in good 
shape; the ablation type gases are somewhat uncertain. 

On Figure 5-53 I am going to touch briefly on the status of 
the experiments. Basically, in terms of laboratory experiments - 
now this is only in terms of the aerothermal environment simulation 
and not the material response - certain aspects can be simulated 
with shock tubes, arcs, lamps, and combined arc-lamp facilities. 
There is no known facility that will do a full job of just cover- 
ing the important parameters that exist. Flight experiment feasi- 
bility studies indicate promise but a lot of expense. It has been 
suggested that we consider shuttle as the launch vehicle which may 
help with the cost problems. 

On the final figure, 5-54, I have selected some priorities 
as to what I think should be done; pretty much in the order that 
I think they should be done, although, for example, the first one 
is certainly just a wish, namely, obtain better input on atmos- 
pheric composition. I am somewhat in agreement with Jerry; I 
think we ought to do some fundamental work in upgrading the tur- 
bulent model. I think we ought to make an effort to continue 
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making the radiation models that are used throughout the industry 
consistent so that we can talk about apples and apples instead of 
apples and oranges. I would like to see something done better with 
the radiation properties for ablation products. I think we ought 
to continue worrying about blowing reduction correlations. That 
is certainly important in terms of planetary entries. I would 
like to see some development of correlations of the inviscid para- 
meters which include the effect of radiation. That capability 
exists; it seems a shame it is not being exploited. I think we 
have to worry the verification tests business further. I would 
like to see some upgrading of the transport property correlations, 
and I think that there ought to be some attention given to the 
non-equilibrium effects. 
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SESSION VI - HEAT PROTECTION 
Chairman: Dr. Phil Nachtsheim 

NASA Ames Research Center 

MR. VOJVODICH: This is kind of like one of the old western 

movies where you can tell the antagonists and the protagonists 
as the guys who wear the black hats and the guys who wear the 
white hats. We have two different view points 'here: the tra- 

ditional approach to the black, carbon phenolic type of heat 
shield and the white, reflecting heat shield. 

DR. NACHTSHEIM: In this session we are going to talk about 

the evaluation of heat shield materials, development of new heat 
shield materials, and then the question of simulation. The 
evaluation will be concerned with the heat shield materials that 
are very well characterized: the carbon phenolic and graphite 

heat shields. Those evaluations will be discussed in terms of 
what was done at the HIP facility in St. Louis and the high- 
powered laser which is here at Ames. In other words, existing 
materials with existing facilities. We will talk about the de- 
velopmental effort on the reflecting heat shield. This concept 
was introduced several years ago, and most people agree that it's 
a good idea. The question remains: how do you do this? So, we 

will be addressing the development of the reflecting heat shield, 
the silica heat shields; and we will have two papers discussing 
that. Then, finally, we will discuss the question of simulation. 
Whether the heat shield be a black heat shield or a white heat 
shield, we do feel that in order to flight qualify it, it should 
be evaluated as closely as possible in the environment that we 
would expect for a planetary entry. 

With that, I would like to introduce Sam Mezines from Mc- 
Donnell-Douglas who will talk about the work he's done on sizing 
the heat shields for Saturn and Jupiter, and some tests he per- 
formed in the HIP facility. 


